Solvable model for electrolytic soap films: 
the two-dimensional two-component plasma 
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We study a toy model for electrolytic soap films, the two-dimensional two-component plasma. 
This model is exactly solvable for a special value of the coulombic coupling constant (3q 2 — 2. This 
allows us to compute the disjoining pressure of a film and to study its stability. We found that the 
Coulomb interaction plays an important role in this stability. Also the adhesivity that measures the 
attraction of soap anions to the boundaries is very important. For large adhesivity the film is stable, 
whereas for small adhesivity a collapse could occur. We also study the density and correlations in 
the film. The charge density near the boundary shows a double layered profile. We show that the 
charge correlations verify a certain number of sum rules. 

PACS numbers: 5.20.Jj, 68.15. +e, 61.20.Qg, 05.70.Np 
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I. INTRODUCTION 



When soap molecules interact with water they disso- 
ciate into anions and cations. The soap anions have 
a hydrophilic head and a hydrophobic tail. Therefore, 
their negative tails prefer to lie on the surface of the film 
while the positive ions can remain throughout the film. 
Soap films have a simple configuration and a well defined 
double layered structure, hence constituting an excellent 
system to be modeled by a Coulomb gas in a confined 
geometry. We are interested in the Coulomb interaction 
between the ions of a two-dimensional soap film and how 
it contributes to the collapse to a black film. 

A soap film, when subject to certain circumstances, 
may collapse to a thickness smaller than visible light wave 
length and therefore it is seen black. This phenomenon 
was observed in the experiments held by O. Belorgey and 
J. J. Benattar |Q and later by D. Sentenac and J. J. Be- 
nattar @. In both cases, salt was added to the soap 
solution. Depending on both the salt concentration and 
temperature, if the external pressure was increased to a 
certain point the soap film collapsed to either a Common 
Black Film (CBF) or to a Newton Black Film (NBF). Ac- 
cording to O. Belorgey and J. J. Benattar M, the CBF 
equilibrium thickness is due to a balance among attrac- 
tive van der Waals forces and repulsive forces between 
the layers. It is believed that short range repulsive forces 
associated to the local structure of water are responsi- 
ble for the NBF stability. The sole difference between 
these two kinds of black films is the thickness of the wa- 
ter layer. The CBF water layer is more than seven times 
as thick as the NBF water layer. But the external lipid 
layer has the same width in both black films. In spite 
of this, the structure and forces involved in black films 
are not completely understood. We want to know if the 
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coulombic forces play a role in this phenomenon. 

In 1997, the mean field Poisson-Boltzmann theory was 
applied by Dean and Sentenac || to three-dimensional 
soap films. Within this framework, they studied the dis- 
joining pressure of the soap film for a wide range of salt 
concentrations and widths of the film. This disjoining 
pressure is the difference between the external and inter- 
nal pressures of the film. However, the phenomenon of 
collapse could not be explained by this mean field ap- 
proach. Later Dean, Horgan and Sentenac Q used a 
functional integral technique to examine a solvable one- 
dimensional Coulomb system model for a soap film. They 
found the film charge distribution and discussed the sta- 
bility criterion for the one-dimensional film. In that 
model they observed the collapse of the film, so one of 
their conclusions was that electrostatic forces play an im- 
portant role in this phenomenon. 

In order to see what aspects are particular to the one- 
dimensional model and what aspects are more general we 
study here another type of solvable model for a Coulomb 
system that can be applied to soap films. Working in 
the framework of classical statistical mechanics, we will 
model the soap film as a symmetric two-dimensional two- 
component plasma, i.e, a neutral system of positive and 
negative particles of opposite charges ±g embedded in 
a neutral background. We are only interested in the 
role that the Coulomb interaction plays in the collapse of 
the film so forces like van der Waals and others will not 
be considered. We will only consider salt-free systems. 
Also, the internal structure of the particles will not be 
regarded. This model is exactly solvable for a tempera- 
ture given by q 2 /HbT = 2. The attraction of the anions 
to the interfaces will be accounted for by a short range 
attractive potential. We will study a film that has an 
infinite surface and a finite thickness. The two dimen- 
sions to be considered lie on the breadth of the film, so 
this system is invariant in one of the two dimensions. We 
will study separately the inner and outer regions of the 
film. For this reason we will use two models. Each one 
is meant to be used to analyze one of the two regions. 

The outline of this paper is as follows. First, in Sec. II 
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we will describe the two models we employ. Then, there 
is brief explanation of the two-component plasma theory 
and of the general method of solution. In Sec. Ill, we use 
the technique presented in Sec. II to find the pressure 
inside the film. Afterward, in Sec. IV, we will show how 



to compute the one-particle densities and the truncated 
two-body densities. Each calculation is followed by its 
corresponding analysis. Finally, we conclude on the role 
that electrostatic forces play in the collapse of this two- 
dimensional soap film. 

The main interest of our work is to study a solvable 
model for a soap film. Given the limitations of this two- 
dimensional model, we can only compare qualitatively 
the structure and behavior of this film to a real one. 



II. THE MODEL AND METHOD OF 
RESOLUTION 

In this section, we will present the two models we em- 
ployed to study a two-dimensional soap film. The two 
models have a lot in common. Both are two dimensional 
systems of particles of charge ±q confined in a slab of im- 
penetrable walls. This aspect is modeled by an infinite 
external potential outside the slab. The anions, nega- 
tively charged, tend to lie on the external surface of the 
film. This is accounted for by an attractive short range 
potential of different form in each model. The cations, 
on the other hand, can lie anywhere in the film and this 
is represented by a constant potential. This potential is 
the same in the two models. 

In the first model (model I) , the short range potential 
is modeled by a delta function, while in model II, it is 
modeled by a step function. The first model will be used 
to find the pressure and the densities in the inner region 
of the film whilst the second model will be employed in 
the computation of the correlations and the densities in 
the outer layers of the film. 

In two dimensions, the Coulomb interaction potential 
of a charge sq at a distance r from another charge s'q is 
logarithmic, of the form v(r) = —ss'q 2 ln(r/d), where d is 
an irrelevant length scale. The adimcnsional coulombic 
coupling constant is T — (3q 2 . For a system of point 
particles the attraction between pairs of opposite sign 
will make the system unstable at low temperatures. For 
this reason the partition function is not well defined for 
r > 2. While for r < 2, the system is stable against 
collapse. 

Let us review the method described by Jancovici and 
Cornu ^| for the two-component plasma. We start with 
the grand partition function 



S= Tn X o / dr i I dr 2 ... ( dr N exp 

N=0 J J J 



(2.1) 



where N is the number of particles and Ao is the con- 
stant fugacity related to the kinetic energy and the 



chemical potential. We consider only neutral config- 
urations: the number of positive particles is equal to 
the number of negative particles. An external poten- 
tial can be described by a position dependent fugacity 
A(r^) = Ao exp (— /3J7 cxt (r,)). In order to avoid diver- 
gences we start with a discretized model. The position 
vector r = (x, y) will be represented by a complex number 
z = x + iy. The particles lie in two interwoven sublattices 
U and V. The positives particles reside in the sublattice 
U with coordinates itj, while the negatively charged par- 
ticles reside in the sublattice V with coordinates Wj. 

For a specific temperature, given by q 2 /ksT = 2, 
this model is exactly solvable. For a configuration with 
N positive particles and iV negative particles, using a 
Cauchy identity, it can be shown that 



exp 



= d 



2JV 



i<j 



det • 



1 



u, 



,N 



(2.2) 



Using this fact, the grand partition function can be writ- 
ten as 



S = det 



1 



d \(vi) 
d A(t) 2 ) d A(ti 2 ) 



d X(ui) d \{u\) 



d AQi) 

El — U2 



Ul—Vl 

d \(u 2 ) 



Ui—V 2 

d \{u 2 ) 



(2.3) 

If each lattice site is characterized by a complex coor- 
dinate z and by a vector which is (1,0) for the positive 
particles and (0, 1) for the negative particles then the 
grand potential can be expressed in the following simpli- 
fied form 



5 = det 



1 



A+(r) \( 
\_(r) {* 







(2.4) 



with A s the fugacity for particles of sign s. 

In the continuum limit where the lattice spacing goes 
to zero (ignoring divergences for the time being) by using 
the identity 
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1 



d 



i 



dz z — z' dz z 
it can be shown that 

-l 



irS (r — r') . 



S = det 



2d z 
2d s 



m + (r) 2d z 
2d z m_(r) 



(2.5) 



(2.6) 



where m s = 2|^A S are rescaled fugacities (S is the area 
of a lattice site). Then defining a new matrix K as 



K = 



2d z 
2d- z 



m + (r) 
m_(r) 



(2.7) 
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the grand partition function H can be expressed as 

5 = det (1 + K) . (2.8) 

The calculation of the pressure, p — —du/dW, reduces 
to finding the eigenvalues of K. While, on the other 
hand, the calculation of the one-particle densities and 
correlations reduces to finding the Green functions G, 
satisfying the following set of equations 

m+(r x ) d xi -id yi \ r( , _ x( ^ 



where 



G ++ 
G_+ G- 



(2.9) 



(2.10) 



and 1 is the unit 2x2 matrix, since it can be shown that 



p Sl (ri) = m sl G sl31 (ri,ri) , 



(2.11a) 

piilT (i"i,r 2 ) = -m sl m S2 G slS2 (ri,r 2 ) G S1S2 (r 2 ,ri) . 

(2.11b) 

When an external potential is acting differently on posi- 
tive and negative particles, it is convenient to define to (r) 
and V (r) as 



m a (r) = m (r) exp [— 2sV (r) 



(2.12) 



To symmetrize the problem for the two types of parti- 
cles, it is useful to define the following modified Green 
functions 

g slS2 (r 1 ,r 2 )=e-^ v ^G slS2 (r 1 ,r 2 )e- s ^\ (2.13) 

Using the following operators 

A = d Xl +id yi +d xl V{r 1 )+id yi V(r 1 ), (2.14a) 
A ] = -d Xl +id yi +d Xl V{v 1 )-id yi V{v 1 ), (2.14b) 

the equations for g ++ and <j__ decouple into 

{to (n) + A 1 " [m (n)] -1 a} g ++ (n, r 2 ) = 6 (r x - r 2 ) , 

(2.15a) 

{771(1-1) + A [m (ri)] -1 A f | g — (ri,r 2 ) = 5 (r x - r 2 ) . 

(2.15b) 

The other Green functions are given by 

ff-+(ri,r 2 ) = - [m(ri)] _1 A 9++ (ri,r 2 ) (2.16a) 
5+ _(n,r 2 ) = [mfnf'AV-fri.rj). (2.16b) 

We will apply the above method to the models explained 
next. 



The film has a thickness W. The outer region has a 
width S, while the inner region has a thickness 2L. The 
breadth of the film is in the x-axis and the film is infinite 
in the y-axis. The origin is set in the middle of the soap 
film. Remembering that m s oc exp (— /3J7 C xt(r)), in model 
I, the position dependent fugacities are 

TO+(r) = to, (2.17a) 
777_ (r) = to + a(S(x - L) + 8{x + L)) . (2.17b) 

The parameter a which we will call adhesivity mesures 
the strength of the attractive potential. 

Model II differs from the preceding in that the attrac- 
tive potential acts over a region of length 5. For a real 
soap film, this thickness is approximately the length of 
the hydrophobic tail. So for this model the position de- 
pendent fugacities are 

m (2.18a) 



TO + (x) 

m_ (x) 



rtii if x G [-L — 6, —L[ U }L, L - 
m if x £ [-L, L] . 



(2.18b) 

with rtii — TOexp(— /3U cx t) > fn since /7 ex t is an attractive 
constant potential. We can therefore distinguish three 
regions: the left border — L — 5 < x < —L (region 1), the 
bulk of the film — L < x < L (region 2) and the right 
border L < x < L + 5 (region 3). 



In this case, the Eqs. (2.15) and (2.16) simplify to 



(m(xi)) - A g±±(x 1 ,x 2 ,l) = m(x 1 )5 (xi - x 2 ) 

(2.19a) 

1 

[~d x T id y ]g±± (x 1 ,x 2 ,l) = g T ± (xi,x 2) l) 

(2.19b) 



m{x\) 

where 
m{x) - 



mo = {rami) 1 / 2 if x in regions 1 or 3 , 
to if x in region 2 . 



(2.20) 

In this model, the potential defined in Eq. ( 2.12| ) is the 
following 

exp(V0 = {(^) 3 |f ^ in regions 1 or 3 , ^ 
I 1 if a; in region 2 . 

With this model we will study the case where mj — > 00 
and 5 — > while keeping their product constant. This 
way model I is a limiting case of model II. 



III. THE PRESSURE 

A. Formal expression for the grand potential 

We shall use here the first model presented in the above 
section (model I), where the position-dependent fugaci- 
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ties are given by Eqs. ( [2.17| ). As explained in the preced- 
ing section the grand potential is given by 



5 = det(l + K) 



(3.1) 



To compute the grand potential we need to find the eigen- 
values of K . The eigenvalue problem for K with eigen- 
values A and eigenvectors (ip, \) reads 



m_(r)x(r) = 2\d E ip(r) , 
m + (r)ip(r) = 2Xd z x(r) . 



(3.2a) 
(3.2b) 



From Eqs. (Q and d2.17|) we find that x is a continuous 
function while ip(x, y) is discontinuous at x — ±L due to 
the Dirac delta distrib ution s in TO _(r). The discontinuity 
of ip is given by Eqs. (3.2a) and (2.17b) 



ip(x = ±L + , y) - ip(x = ±L ,y) = -\{x = ±L, y) . 



(3.3) 



Inside the film, for — L < x < L, Eqs. (3.2) can be com- 
bined into the Laplacian eigenvalue problem 



/ m \ 2 

a *=(a) 



(3.4) 



Due to the translational invariance in the y-direction we 
look for solutions of the form 

X(r) = (Ae- K * x + Be K ' x )e ikv , (3.5) 
where k* = (fc 2 + (to/A) 2 ) 1 / 2 . From Eq. ( |3.2b| ) we find 



fl>(r) 



X 



(A(k - K*)e~ K ' x + B{k* + k)e K " x )e 



hy 



Outside the film Eqs. (3^2) reduce to 

d s ip = and d z x = . 



(3.6) 



(3.7) 



That is ip is analytic and \ is antianalytic. Since we are 
looking for solutions with y dependence e lky this gives 



Ce kz = Ce kx+tky , 

£) g -fcz _ jj e -kx+iky 



(3.8a) 
(3.8b) 



In order to have vanishing solutions at infinity, from the 
preceding equations it is necessary that for k > 



and for fc < 



V>(r) 
X(r) 



V-(r) 
X(r) 



for x > L , 
for x < —L . 



for x < —L . 
for x > L . 



(3.9a) 
(3.9b) 



(3.9c) 
(3.9d) 



Eqs. (3.3) and (|3.9|) are the boundary con ditio ns that 



complete the Laplacian eigenvalue problem (3.4). These 
boundary conditions yield a homogeneous linear system 



for the coefficients A and B, which in the case k > 
reads 



A B 



(A 2 (fc — K* 

(\ 2 (n* + k)+am)e K 



am)e K ' L e K ' L 



(3.10) 



In order to have non trivial solutions the determinant of 
this linear system must vanish. This gives the following 
equation that must be satisfied by the eigenvalues A 



{\ a {k* + fc) + am)e 2K * L + (A^(k* - fc) - am)e 



-2k* L 



0, 

(3.11) 



that can also be written as 



cosh(2K*L) 



sinh(2K*L)/ K * = 0. (3.12) 



A similar equation is found for the case fc < in which 
one should change fc for — fc. As a consequence of this fact 
the set of solutions for fc < is the same as for fc > 0. 
From now on we will only consider the case fc > 0. 
The grand potential per unit length is then given by 



(3w = 



1 f 

-7T / lnJJ(l + A fc )dfc 

1 f +OC 

— / inTT(i + A fe )dfc, 

n Jo 



(3.13) 



where the product runs over all A& solution of Eq. (3.12). 
This product can actually be performed as explained in 
refs. |6|, 0, B|. Let us introduce the analytic function 



/*(*) 



cos h(2v/fc 2 +m 2 z 2 L) 

sinh(2Vfc 2 + m 2 z 2 L) 



+ (fc + amz 2 ) 



\Jk 2 + m 2 2; 2 



(3.14) 



By construction the zeros of fk are the inverse of the 
eigenvalues A^. This function fk can be factorized as a 
Weierstrass product running over its zeros. Since / z (0) = 
exp(2fcL), /'(0) = 0, f(z) = f(-z) and the zeros of f k 
are 1/Afc, the Weierstrass product representation reduces 
to 



f k (z) = exp{2kL)]J{l - z\ k ) 



(3.15) 



Then the product appearing in the grand potential ( 3.13j ) 
is simply /&(— l)e~ 2kL . Finally, the grand potential per 
unit length reads 



dk 



In I cosh(2Ki) 



- 2fcL + 
fc4 



(3.16) 



■ sinh(2Ki' 



where k = (fc 2 + to 2 ) 1 / 2 . The above integral is actually 



divergent and should be cutoff to a fc n 



1/R where R 



is the diameter of the particles as expla ined in Ref . || . 
It can be checked that for a = 0, Eq. ( |3.16 ) yields the 
known grand potential for a two-component plasma in a 
strip of hard walls ((| . 

It is interesting to look at the large- L behavior of the 
grand potential, 



-W Pb + 2 7 + O (e- mW ) 



(3.17) 



where W = 2L is the width of the film, the bulk pressure 
is given by 



1 



7T 



l/R 



(k — k) dk 



i) 



2ir V n mR 



(3.18) 



(a known result from Refs. |5|, in the limit of vanishing 
cutoff R — > 0), and the surface grand potential is 



i 

2^ 



l/R 



In 



mot 



dk. (3.19) 



In the limit R — > the surface grand potential reads 



07 



in 

47T 



a In 



mi? 



1-7T- 



1-a 2 



a - 



ln(a + i; 



(3.20) 



When a = the above expression reduces to the known 
result H 1 



/3 7 (" = 0) 



m 
27 



1 



(3.21) 



When the cutoff R vanishes the surface grand potential 
diverges (except for a = 0). This is expected since neg- 
ative particles are strongly attracted to the boundaries 
and for point particles this would create divergences in 
the surface grand potential in addition to the usual diver- 
gences in the bulk pressure due to the collapse of particles 
of opposite sign. 



Finally it should be noted in Eq. (3.17) that there are 
no algebraic corrections in 1/W to the grand potential. 
The next correction after the surface term is exponen- 
tially small. This is the same situation as in a strip with 
hard walls (a = 0) but very different from the situa- 
tion of ideal conducting boundaries Q and ideal dielec- 
tric boundaries Jl(|. In those later cases there is indeed 
an algebraic universal finite-size correction to the grand 
potential equal to 7r/24VF". 



B. The disjoining pressure 

The pressure in the film can be obtained from the 
grand potential as p = —duj/dW. The disjoining pres- 
sure is defined as the difference between the pressure of 
the film and the pressure of an infinite system (W — > oo, 




FIG. 1: The disjoining pressure pd as a function of the 
width W for several values of a. From top to bottom 
a = 2,1,0.5,0.3,0. Notice that the disjoining pressure be- 
comes positive for small W except in the case a = where it 
is always negative. 



the bulk pressure) : pd = p — Pb_M ■ Using Eq. ( 3.16| ) for 



the grand potential and Eq. (3Tq) for the bulk pressure, 
we find 



(3pd = - 



g(k) dk, 



with 



g(k) = 



2k(—k + k + ma)e 



-2kW 



K + k + ma + (k— k — ma)e 



-2kW 



In the limit W — > 0, the total pressure is 



/3p(W 



ma 



(3.22a) 



(3.22b) 



(3.23) 



For W = and a^0, the pressure diverges as l/R when 
the cutoff vanishes, stronger than the usual logarithmic 
divergence. From this fact, it is clear that for a ^ the 
disjoining pressure will be positive when W — > (actually 
Pd — > +oo as W — * 0) . The case a — is particular since 
then the disjoining pressure is of same order as minus 
the bulk pressure for small- IF and then pd — > — oo when 
W -> 0. 

For non-zero width films W ^ the disjoining pressure 
is finite for vanishing cutoff R — * (this limit has already 
been taken in Eq. (3.22)). 

Fig. (|]) shows several plots of the disjoining pressure as 
a function of the width W, for different values of a. We 
notice two different behaviors. For a = 1 and a = 2, the 
pressure is a monotonous decaying function of the width. 
This indicates that the film is stable for all widths. For 
a = 0.3 and a = 0.5, we notice that the pressure is no 
longer a monotonous function of the width. There exists 
a "critical" width W c and for W > W c the pressure is 
an increasing function of the width. This indicates that 
the film is unstable for W > W c . Indeed, in that region 
a small change in the applied pressure to the film will 
result into a collapse to a width smaller than the critical 
W c . This transition is discontinuous (first order). This 
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FIG. 2: The disjoining pressure pd as a function of the width 
W for a = 0.3. The critical length is W c . The phenomenon of 
collapse is illustrated as follows. Initially the film has a width 
Wa corresponding to the point A in the pd - W diagram. A 
small change in the pressure will force the system to go to 
point B: the film has collapsed to a width Wb < W c < Wa- 
This transition is clearly first order (discontinuous). 




FIG. 3: The critical width W c of the films as a function of 
the adhesivity a. Films larger than W c are unstable. 



is illustrated in Fig. (||). In the case a = the disjoining 
pressure is always negative and increases if W increases. 
The film is unstable for all widths when a = 0. 

The critical value of a distinguishing between these two 
different behaviors is a c = 1 as it will be shown below. 
In order to determine if the pressure is an increasing or 
decreasing function of W we study the sign of dpa/dW 
and in particular the sign of dg(k) / dW . We have 



2/; IV 



dg(k) mK 2 [m(l - a 2 ) - 2ak] e 

dW [k + k + ma + (k - k- ma)e- 2KW ] 



(3.24) 



Clearly for a > 1 the function dg(k)/dW is negative 
for all values of k > 0, and consequently the disjoining 
pressure will be a decreasing function of W. 

For a < 1 the function dg(k)/dW is positive for values 
of k < k* — m(l — a 2 ) /2a and negative for values of k > 
k*. Since for large values of W, the function dg(k)/dW 
decays exponentially, the dominant part of the integral 



will collapse, mimicking the collapse to a Common Black 
Film or to a Newton Black Film in real soap films. 

This situation is somehow different to the one exposed 
in Ref. [Q for a one-dimensional film. Common features 
of the present study and the one-dimensional case are 
that for sufficient large values of the adhesivity a a sta- 
ble film region exists. For small values of a (a < 1 in 
our case) a collapse can occur. But the main difference 
is that for a < 1 very large two-dimensional films are un- 
stable while very large one-dimensional films are always 
stable for a > 0. Another important difference is that in 
the one-dimensional case multiple collapses are possible 
whereas in the two-dimensional case we only have one 
collapse (or no collapse). 



d(3pd 
dW 



dg(k) 
dW 



dk . 



(3.25) 



will be given by small values of k, where dg(k)/dW is 
positive. Then, for large values of W, dpd/dW will be 
positive and pd will be an increasing function of W for 
large W. 

The exact value of W c where dpd/dW changes of sign 
cannot be determined analytically in a simple manner. 
However, it can be determined numerically. In Fig. (|^) 
we plot W c as a function of the parameter a. 

As a general conclusion of this analysis it can be said 
that the attractive potential in the boundary of the film, 
whose strength is characterized by a, allows the film to 
stabilize. For a > 1, films of arbitrary width are stable, 
whereas for a < 1 only small films are stable, larger films 



IV. DENSITY AND CORRELATIONS 



The density and correlations can be obtained by com- 
puting the Green functions introduced in Sec. II. The 
present section is divided into two parts. First, we will 
study the density and correlations inside the film using 
model I. In the second part, we will compute the density 
in the boundary of the film and the correlations when 
one point is on the boundary of the film. For that part 
we will need to use model II since the Green functions of 
model I are discontinuous on the boundary as we will see 
below and therefore do not give any information about 
the boundary. 
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Inside the film 



1. The Green functions 



In the present geometry it is natural to work with the 
Fourier transform G ss i of G ss > in the y-direction 



Gw(ri,r 2 ) 



dk 

G ss ,( Xl ,x 2 ,k)e ik ^-v*^. (4.1) 
Zir 



Then Eq. (2.9) translates in Fourier space to 
m+(xi) d Xl + k 



d Xl - k m-(xi) 



G(xi,x 2 ,k) = S(xi - x 2 )l ■ 



(4.2) 



Let us detail the calculation of G__ and G-\ , the one 

for G++ and G |_ follows similar steps. The equations 

are 

toG+_ + (d Xl + k)G— = 0, (4.3a) 
(<9 Xl - k)G+- + m_(zi)G__ = 5(xi - a; 2 ).(4.3b) 



From Eq. (4.3a) we deduce that G is continuous. How- 
ever because of the Dirac delta distributions in the defi- 
nition of m_ the function Gj will be discontinuous at 



X\ = ±L. From Eq. (4.3b) we deduce the discontinuity 
of G^ at x\ = ±L if X2 =f= x\ 

G+_(xi = ±ZT) - G + -{ Xl = ±L + ) = aG—{xi = ±L) , 

(4.4) 

If both points i"i and r 2 are inside the film but not on 
the boundary both fugacities are equal m + = rn_ = m 



and then Eqs. (4.3) can be combined into 

{0 2 X1 - k 2 )G—(x 1 ,x 2 ) = -mS( Xl - x 2 ) , (4.5) 



with k = (k 2 + m 2 ) 1 / 2 , while G^ is given by Eq. (4.3a). 

If ri is outside the film while r 2 is fixed inside the film , 
thenm+(xi) = m-(x\) — and the solution of Eqs. (4.3) 
is 



G + -(x u x 2 ,k) = De kxi . 



(4.6a) 
(4.6b) 



In order to have finite solutions at X \ = ±oo it is neces- 
sary that for k > 

G (xi < —L) = and G+_(xi > L) = ,(4.7a) 

and for k < 

G__ (xi >L) = and G+_ (x x < —L) = .(4.7b) 



Eqs. (4.4) and (4.7) are the boundary conditions that 
complement the differential equations ([L3J) for the Green 
functions. 

The solution is of the form 



G — (xi,x 2 ) = £![ e -»l«*-*»l + 



m 
2k 1 

AeT KXl + Be KX2 ] 



(4.E 



where the coefficients A and B are determined by the 

boundary conditions. Finally the Green function G 

is, for k > 0, 



J 



G — (xx,x 2 ,k) 



2 k 



g — k\xi— Xl\ _j_ 

-(k + k + am)e~' i< ' Xl+x ^ + (k - k - am)e K{xi+x ^ 



+ 



(n-k- am)e- KW + (k + k + am)e KW 

2(—n + k + am)e~ KW cosh k(xi 



x 2 ) 



(k — k — am)e kW + (k + K + am)e 



kW 



(4.9a) 



and for k < 0, 



G — (xi,x 2 ,k) = — 



-K.\X1 —X 2 



(n + k- am)e-*( I1+ls ) + (—k + k — am)e K( - Xl+x ^ 
(n + k- am)e~ KW + (« - k + am)e KW 

2(—k — k + amie 



+ 



cosh k(xi — x 2 ) 



(k + k — am)e kW + (n — k + am)e K 



w 



(4.9b) 
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and G^ can be obtained from Eq. (4.3a). Similar calculations lead to G++ for fc > 0, 



G ++ (x 1 ,x 2 ,k) = — 
2k 



while for fc < 0, 



+ 



-(re + k) (1 - a^) e ^ Xl+X2 '> + (re - k) (l + a^) e^ 1 ^ 
(re — fc — am)e _KVI/ + (« + fc + ar7i)e KH/ 

2(fc — k + am)e _KW/ cosh k(x\ — X2) 



+ 



(k - k - am)e- KW + (re + fc + am)e* 



(4.10a) 



G ++ (xi,a;2,A:) 



2re 



e ~K\xi-X 2 \ , 



(re + fc) (l + a^) e^* 14 - 2 ^ — (re — fc) (l — e -"(*i+*») 
(ft + k - am)e~ KW + (re - k + am)e KW 

2(—k — k + am)e~ KW cosh k(xi — X2) 



(« + k — am)e kW + (ft — k + am)e K 



r 



(4.10b) 



The Green function G |_ is obtained from 

G- + (xi,x 2 ) = ~—(d xl - k)G ++ (xx,x 2 ) ■ (4.11) 
m 

We omit the details of the calculation for G++. Let us 

only note that in this case, G^ is continuous while G++ 

is discontinuous at x\ = ±L with 

G ++ { Xl = ±L~) - G++O1 = ±L+) = aG+_(xi = ±L) . 

(4.12) 

The Green functions G(i"i,r2) in po sitio n space are 
given by the Fourier transform formula (4.1). The term 
mexp(— ft|xi — X2|)/(2ft) that appears in all G ss give the 
bulk contribution to G ss H 



Gbuik = —K (m\r 1 - r 2 |) 
Ztt 



(4.13) 



where Kg is the modified Bessel function of the second 
kind of order 0. 

For a — the expressions for the Green functions re- 
duce to known results |ll|. 



2. The density 



The density of species of charge s is given by 



p s (r) = m s (r)G ss (r,r) 



(4.14) 



Using Eqs. ( |4.9| ) and ( 4.10|) for the Green functions we 
obtain the densities 



z (fc — ft + am)e kW — (k + am) cosh( 2 kx) dk 

p-(x) = Pb + — - ( r \ _ kW \ TTX hw 1 — ( 415a ) 

7T Jo (ft — fc — am)e KVV + [k + k + am)e KVV 



m 2 poo (fc — ft + am)e kW — ( fc — am — ^Jl. j cosh(2ft.x) ^ 
7r J a (ft — fc — am)r* + (ft + fc + am)e K 



P+ (x) = Pb + — j — „_^_ kW ,,.,,s kW — (4.15b) 



where P b is the bulk density (actually divergent when (measured in units of q) is then 
the cutoff R — s- 0). The charge density p = p + — p_ 

, , 2ma f°° Ke~ KW cosh(2rex) dk 

p(x) = — 



k + fc + am + (k — fc — am)e 2kW 

(4.16) 
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FIG. 4: The charge density as a function of the position x 
for several values of a. The width of the film is W — 8/m. 
From top to bottom a — 2,1, 0.5, 0.3. 



Fig. (|J) shows several plots of the charge density as a 
function of the position x. This figure can be understood 
as follows. Because of the strong attractive potential on 
the boundary an important part of the negative particles 
(the soap molecules) are stuck in the borders of the film 
at x = ±L creating a layer of negative surface charge den- 
sity (actually it is really a linear charge density since our 
system is two-dimensional, in a three-dimensional case it 
would be a real surface charge density) . In the framework 
of model I, this negative surface charge density cannot be 
seen in Fig. (Q) nor in the analytic expressions found for 
the densities, but it will be stu died i n detail in the sec- 
ond part of this section (section IV B| ) when we will work 
model II. 

The system on the interior of the film is then non- 
neutral with an excess positive charge. This excess pos- 
itive charge screens the negative surface charge density 
at the borders as it can be seen on Fig. (|]). The den- 
sity near the boundary is positive and becomes very large 
when the boundary is approached. Away from the bor- 
ders, near the middle of the film, the system is almost 
neutr al. A lso in Fig. (||) and from the analytic expres- 
sion (4.16) for the density it can be checked that the 
screening length is of order m _1 , a well-known result [^J. 
It should be noted that for a — 0.3 and a = 0.5 the 
system should collapse according to the analysis of the 
preceding section (Sec. III). However there is no hint on 
the charge density profiles indicating the collapse. This 



was also the case on the one-dimensional model M . 

The surface charge on one border —a can be computed 
by using the screening sum rule 



p(x) dx , 



(4.17) 



giving 



2tt 



(1-e 



-2kW\ 



dk 



K + k + am + (k — k 



am)e~ 2KW 

(4.18) 



Actually the above expression is divergent and should 
be cutoff to a fc max ~ 1/ R as it has been done for the 
pressure. In section IV B a more direct calculation of a 



will be done and the validity of the screening sum rule 
will be proven. 

The surface charge density a is an increasing function 
of W. Actually for large films it converges exponentially 
fast to the value 



am 
m 

4-7T 



1/R 



1 





a In 



k + k + am 
2 a 2 + l 
mR 



a 



dk 
hi(a 



D + i 



.(4.19) 



The vanishing terms when R — > have been omitted in 
the second equality. The dominant term of <7b, when the 
cutoff vanishes, is proportional to a. It is clear that a 
controls how much the boundaries get charged. 

For very large films W — > oo, it is interesting to study 
the relationship between the surface charge density Ub 
and the surface tension. Actually the surface grand po- 
tential 7 computed in Sec. Ill is really the surface tension 
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of the system since we are dealing with polarizable inter- 
This is clear in model II (model I being 



faces 

a special limit of model II) where the particles are free 
to go from region — L < x < L to the boundary re- 
gions — L — 5 < x < —L and L < x < L + S. In that 
model the control parameter for charging the boundaries 
is the attractive potential J7 GX t , or equivalently the fugac- 
ity m.j. For model I the control parameter is the adhe- 
sivity a. From the forma l expressions of 7 and (T{, given 
by Eqs. fl3.19| ) and ( [4.19D it can be verified that 



cr h = -pa— , 
oa 



(4.20) 



which can be regarded as a particular form of Lippmann 
equation for this model 



3. The correlations 

The truncated two-body correlation function for a par- 
ticle of sign s at 17 and a particle of sign s' at r2 is given 
in terms of the Green functions by 



Pss> ( r l' r 2) 



-m s (ri )m a > (r 2 )G ss > (17 , r 2 ) G s > s (r 2 , 17) . 

(4.21) 



With the expres sions for the Green functions given by 
Eqs. ( |4.g| ) and (4.10) the correlation functions can be 
obtained. 

Due to the screening properties the structure near a 
boundary will not be modified considerably by the pres- 
ence of the other boundary. For this reason it is inter- 
esting to study in further detail the case of large films 
W — > 00. The corrections for finite films are exponen- 
tially small in mW. 
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For W — > oo, taking the origin at the boundary (x 
is now changed to x + L), the Fourier transforms of the 
Green functions simplify, for k > 0, to 



G — (x 1 ,x 2 ) = G hu i k -—e 

2k 



— k(xi+X2) 



(4.22a) 
(4.22b) 



G ++ (xi,x 2 ) = Gbuik + 

, m (tt-fc) [1 + £(*-*)] p - K( , 1+ , 2 ) 



and for fc < 

G — (3:1,352) = 

6+4.(3:1,352) : 



2k k + k + am 



Gbulk + 

to ac + fc - am e _ K(xi+a;2) 
2ac ac — fc + am 

Gbulk + 

TO (fc Av) [!_£(«_*)] 



(4.22d) 

-«(a;i+X2) 



2k 



At — fc + am 



As a very curious fact it should be noted that the pre- 
ceding expressions for G are formally identical with 

the ones of Ref. || for a two-component plasma near 
a charged hard wall if one chooses in the later case 
the external surface charge density of the wall equal to 
—am/(2ir). This does not mean that in our case the 
charge density (which is actually not external to the sys- 
tem, but internal and due to the reorganization of charges 
in the him) is —am/(2ir). In the preceding subsection 
we have computed the surface charge density —a and we 
know that it is not equal to — am/ (2tt). Furthermore for 
the Green functions G++ this comparison does not hold, 
the expression for G++ given by Eqs. ( 4.22 ) and those 
from Ref. § are very di fferen t. 



It is clear from Eqs. (4.22) that the correlation func- 
tions will decay exponentially fast in the n-direction 
(through the width of the film) . However it is well known 
that in the y-direction, parallel to the boundary, the cor- 
relation functions usually decay algebraically [|l3| [l4|, [l5| . 
For two-dimensional Coulomb systems, near a hard or 
dielectric (non-conducting) wall, they should decay as 
y~ 2 . For a hard wall (a = 0) the total charge correla- 



(2)T 
P-+ 
V2 -> 



(ri,r 2 ) 




tion function S{r 1 ,r 2 ) = /3+ + T (ri,r 2 ) + (ri, r 2 ) 



ri,r 2 ) should behave, for y = y% 



S(ri,r 2 ) 



f{xi 1 x 2 ) 



(4.23) 



and the function f(x±, x 2 ) should obey the sum rule 



I dX2 l dx 1 f(x 1 ,x 2 ) = -- 2p 



(4.24) 



The preceding asymptotic behavior is very general for a 
Coulomb system near a plane hard wall. 

Returning to our case, the Fourier transform of the 
Green functions have a discontinuity at fc = which 
will translate in position space into a decay as 1/y for 



large y, then the correlation functions will indeed have 
a decay as 1/y 2 . More precisely, the large- y behav- 
ior of the charge correlation function is found to be 
S(ri,r 2 ) ~ f(x\,x 2 )/y 2 with the function f(x\, x 2 ) given 
by 



f(xi,x 2 ) = -- 



,2.-2111(11+12) 



(a + l) 2 7T 2 
This function obeys the sum rule 



poo /* 00 

(4.22c) J dx z J dx 1 f(xi,x 2 



1 



27r 2 /3(a + l) 2 



(4.25) 



(4.26) 



According to Ref. |13|, [14| |15j] for a Coulomb system near 
a plane hard wall w ith an external charge density on the 
wall the sum rule ( [1.24 ) is not modified. In our case, 
where the boundary is charged by a fraction of the par- 
ticles of the syst em, the sum rule is modified. However 



in the sum rule ( 4.26 ) we only accounted for the correla- 
tions of particles in the fluid. As we will see in Sec. IV B 



there are also correlations for particles that are absorbed 
in the bound ary an d when these are taken into account 
the sum rule ( 4.24 ) is verified. 

It is also interesting to comment on the case a — > 00. 
It can be checked that in that limit Eqs. (4.22) reduce to 



G — (xi,x 2 ) = Gbuik-r-e 

Zk 



G ++ (xx,x 2 ) 



Gbulk + 
, (*~fc) 2 

2mK 



(4.27a) 
(4.27b) 



-k(xi+x 2 ) 



for all values of fc. Computing the inverse Fourier trans- 
forms gives 



G__(r 1; r 2 ) = — [K (mr 12 ) - K (mr* 12 )](4.27c) 



G++(ri,r 2 ) = —K (r 12 ) + 



(4.27d) 



rn 



+ ^e-^K 2 (mr* 12 ) 1 



where r\ 2 = |ri — r 2 | and r* 2 = |ri — r|| with r 2 — 
(—x 2 ,y 2 ) being the image of r 2 . The angle 4>* 2 is the 
angle of the vector ri — r' 2 with respect to the ir-axis. 

It is clear that, when a — > 00, the Green functions have 
no longer an algebraic decay along the y-direction. The 
decay is now exponential in all directions and of the form 
cxp(— TOJ-12) and exp(— mr* 2 ). 



B. On the boundary of the film 

We are now interested in the structure of the film at the 
boundary. Model I cannot give directly any information 
on the density or the correlations at the boundary since 
some of the Green functions are discontinuous there. We 
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shall use instead model II where the fugacities are given 

by 



m-(x) 



m l ifxe[-L-S,-L[U]L,L + S], 



if x e [-L, L] 



(4.28) 



and m + (x) = to everywhere. Recall that we distinguish 
between three different regions: the left border — L — 5< 
x < —L (region 1), the bulk of the film — L < x < L 
(region 2) and the right border L < x < L + 5 (region 3). 



1. The Green functions 

It is useful to work with the modified Green functions 
as explained in Sec. II. As before we will concentrate on 

the computation of g , the one for g ++ follows the same 

lines. We fix the source point r 2 in region 1. Then, the 
Green function obeys Helmoltz equation in the different 
regions 



(A ri - m{x 1 ) 2 )g — (ri, r 2 ) = -m(x 1 )d(r 1 - r 2 ) . 



(4.29) 



and 



ff+-(ri,r 2 ) = -. , {-d Xl +id yi )g—(r 1 ,T 2 ), (4.30) 
with the position dependent fugacity 



m(x) 



to 



= {mmi) 1 / 2 if x in regions 1 or 3 , 



if x in region 2 . 



(4.31) 



Working with the Fourier transforms g ss i of g ss i we have, 
if ri is in region 1, 



mo 

+A 1 e- K « Xl + B ie KoXl 



g—(x 1 ,x 2 ) = -° e-" ^ 1 -*^ + (4.32a) 

K 



if ri is in region 2, 

g—{xx,x 2 ) = A 2 e~ KXl + B 2 e KXl , (4.32b) 

and if ri is in region 3, 

g—{x u x 2 ) = A 3 e- K ^ + B 3 e K ^ , (4.32c) 

with Kg — (toq + k 2 ) 1 / 2 . The coefficients Ai and Bi are 
determined by the following boundary conditions: G ss i 
should be continuous at x\ — ±L and at x\ — ±(L + S). 
Furthermore, for k > 0, G__ = if x\ < —L — 5 and 
G+_ = if x x > L + S. And for k < 0, G__ = 
if x x > L + 5 and G+_ = if xi < -L - S. These 
boundary conditions yield a linear system of six equations 
for the coefficients Ai and Bi for each case k > and 



k < 0. This system can be solved by standard matrix 
manipulation programs like Mathematica. Since the 
solution is actually not very illuminating and too long 
to reproduce here we will only consider from now on the 
limit Tfii — > oo, 5 — > while niiS = a is kept constant. 
This limit is taken after the linear system has been solved. 
In this limit model II reduces to model I. 

In the limiting procedure it is very important to ob- 
serve how the Green functions scale with the fugacity m^. 

1 12 

We find that is proportional to to/ when ri is in 
region 1. Then the density will be proportional to rrii 
and diverges in the limit rrii — > oo. However the "sur- 
face" charge density er_ = 5p~- will have a finite value. 
Similar scaling behaviors appear in the other regions, giv- 
ing finite surface charge density-surface charge density 
(cr_(ri)cr_(cr(r 2 )) correlations for both points in bound- 
ary region 1 , and finite charge density-surface charge den- 
sity (p s (ri)cr_(r 2 )) for ri in region 2 and r 2 in region 1, 
as it should be. 

On the other hand the function g ++ for instance is of 

1/2 

order l/m i when ri and r 2 are in region 1. This will 
give a finite charge density p + in the boundary which is 
negligible when compared to the surface charge density 
cr_ . The same holds for correlations when ri is in region 
2 and r 2 in region 1, we find finite charge density-charge 
density correlations (p s (ri)p + (r 2 )) which are negligible 
when compared to the charge density-surface charge den- 
sity correlation function (p s (ri)<r_(r 2 )). 

This is not surprising. We know from the preceding 
section that the strong attractive potential has created a 
negative surface charge density at the boundaries. 

The results for the relevant Green functions are, for ri 
in region 1 and k < (with W = 2L) 



9— = 



m {l~e- 2KW ) 



k — k + am + (k + k — am)e 



-2kW 



, (4.33) 



while for k > 0, = 0(to,q 1 ). All other Green func- 
tions in this region are of order 0(m,Q 1 ). 
For ri in region 2 we find, for k < 0, 



9—( x i) = 



2(m m) 1 / 2 e- KW sinh k(L - Xx) 
k — k — am + (k + k — am)e~ 2KW 

(4.34a) 



and 



9+-{x\) 



to 



1/2 



(4.34b) 



-k(l+x 1 ) [ K -k + {n + k)e~ 2 < L - x ^] 
k — k + am + (k + k — am)e~ 2KW 



while for k > the Green functions g__ and g^ arc 

x/ 2 

of order O(m ). Also the other Green functions g++ 
and g — I- are of order 0(to,q 1 ^ 2 ). 

When rj is in the other boundary (region 3) we found 
that all Green functions are of order 0(to,q 1 ). 
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2. The density 

W e now compute the density on the boundary. Using 
Eq. (4.33) for the Green function g and the formula 



p s {r) = m g ss (r,r) , 



(4.35) 



we find a charge density p_ proportional to rrij. Then 
the surface charge density cr_ = bp- is finite in the limit 
mi — > co, b — > and mi<5 = a fixed. We have 



2- 



(1 



(/A- 



/c — k + am + (k + k — am)e 2kW 

(4.36) 



Making a change of variable k — > — fc in the integral 
we find that the surface charg e density cr_ is equal to 
the one computed in Sec. IVA using the screening sum 
rule ( 4.17 ), <7_ = a, as it should be. The screening sum 
rule (4.17) is then verified. 



3. The correlations 

For both points ri and r 2 on the boundary (region 1) 
we compute the charge density correlation function by 
using 



p\'_ (ri,r 2 ) 



-ml\g— (rr,r 2 )| 2 



(4.37) 



The Green function g given by the Fourier transform 

of Eq. ( |4.33] ) is proportional to mo = ^Jmmi. Then, the 
charge density correlation function is proportional to m 2 . 
This gives a well defined surface charge density-surface 
charge density correlation function 



(a-( Vl )a-(y 2 )f = b 2 p ( ? ) I(r 1 ,r 2 ), 



(4.38) 



in the limit m, — > 00, b — > and niib = a fixed. The 
final expression for the correlation function is 



(cr_(yi)f7_(y 2 )) J 



ma 



X 



(1 



. 2tt 

2KWyky dk 



(4.39) 



k + am + (k — k — ma)e 2kW 



with y = yi — y 2 . As before it is interesting to study the 
decay of the correlations along the y-axis. The disconti- 
nuity of the Fourier transform of the Green function at 
k = will be translated into a 1/y decay. Then the sur- 
face charge correlation will have the asymptotic behavior 



(cr_(yi)cr_(y 2 )) J 



(1 



-2mW\2 



1 



4?r 2 (1 + a + (1 - a)e- 2mW Y V 1 



and for very large films W 
(cr_(?/i)cr_(2/ 2 )) T ~ - 



00 



1 



47r 2 (a + l) 2 y 2 



(4.40) 



(4.41) 



In relation to the results of Sec. IV A on the asymptotic 
behavior of the total charge correlation function S'(r 1 , r 2 ) 
inside the film we notice that for y — > 00 and in the limit 
W — >■ 00 the following sum rule is verified 

roc roc 

(°"-(yi)c-(y2)) T = a 2 / dx 2 dx 1 S(r 1 ,r 2 ) ■ 
Jo Jo 

(4.42) 

When ri is in region 2 (inside the film) and r 2 remains 

(2)T 

on the boundary, the correlation function p__ is given 
by 



P— (ri,r 2 ) 



-m m| 5 __(ri,r 2 )| 2 , (4.43) 
with the Green function g given by the inverse Fourier 



transform of Eq. ( |4.34a| ). The correlation function /9+l T 
is given by 



p+l T (ri,r 2 ) = m TO|5+-(ri,r 2 )| :i; , 



(4.44) 



with ff-| given by the inverse Fourier transform of 



Eq. ( 4.34b| ). We notice that both correlation functions 



are proportional to m 2 , = mim. Then the charge density- 
surface charge density correlation (/3 s (ri)cr_(y 2 )) = 
( 2 \ r p 

5p s _ (ri,r 2 ) will have a well defined finite limit when 



m, — > 00, S - 
\\ r find 



and rriiS = a. 



(p_(ri)o-_(y 2 )) = 



(4.45) 



e~ KW sinh [k(L - xi)] e iky dk 



k + k + am + (k — k — am)e 



-2kW 



and 



(p + (r 1 )a-(y 2 )) J 



47T 2 



k + k + am + (k — k — am)e~ 2liW 



(4.46) 



dk 



It is interesting to notice a relation between 
(p- (i"i)<t_ (y 2 )) T and the surface charge density 
correlation when both points are at the border, which 
we might call "continuity" . This relation is 

(p-{x! = -L, yi )a^(y 2 )) T = — (o-_(yi)CT_(y 2 )) T . 

a 

(4.47) 

It is also important to study the decay of the correla- 
tions along the boundary, when \y\ = \yi — y 2 \ — * 00. 
Here again the discontinuity of the Fourier transform of 
the Green functions at k — is responsible for a 1/ y 2 de- 
cay of the correlations. The asymptotic behavior of the 
correlations for |yj — y 2 \ — > 00 is 



(P-(ri)cr^(y 2 )) J 



am [sinh m(L — x\)] 2 e 2mW 

27, 2 ' 



[1 + a + (1 - a)e- 2mW ] n V 



(4.48) 



13 



and 



(p+Oi)o--(2/2)) J 



am [cosh m(L — xi)] e 2mW 

[1 + a + (1 - a)e- 2mW f *V 

(4.49) 



contains the contribution when both particles are on the 
boundary. The origin is taken on the boundary. Using 
the asymptoti c exp re ssion s for the d ifferent correlations 
given by Eqs. ( 4.25 ), ( |4.5C| ) and ( 4.41 ) one can check that 
the sum rule (4.24) is verified: 



For very large films W — ► oo, taking now the origin at 
the boundary x\ — ► x\ + L, the total structure function 
(p<r) T = (p_(T_) T — (p + <j-) T has the asymptotic behav- 



dxt / da;2 5 t otai(ri,r 2 ) 
ii Jo 



(p{r 1 )a{y 2 )) J 



—2mxi 



2n 2 (a + l) 2 y 2 



(4.50) 



1 1 + 2a + a 2 
47r 2 y 2 (a + l) 2 

= - w ■ (455) 

For finite W there exists also a sum rule similar to 



We notice that this asymptotic correlation function obeys 
a sum rule with the charge correlation function S(ri, r 2 ) 
for both points inside the fluid 



Eq. (4.24) for conducting Coulomb systems 
sum rule reads 



15 



/L pL 
dx 2 / cfeiS , totai(ri,r 2 ) 
-L J -L 



1 



I3ir 2 y 2 



r§. This 



(4.56) 



S(ri, r 2 ) dx 2 = a (p{v 1 )a(y 2 )) J 



(4.51) 



and also a sum rule with the surface charge density cor- 
relation when both points are on the boundary 



for \y\ — > oo. From Eqs. ( |4.9| ) and (4.10) one can com- 
pute the asymptotic behavior of the correlations for both 
points in the fluid. We find 



(p(ri)cr(y 2 )) dx x = a {(r(y 1 )a{y 2 )) J 



(4.52) 



5(ri,ra) 



cosh [2m(xi + x 2 )\ 



2n 2 y 2 [cosh(mFF) + asmh(mW)Y 



The remaining case to complete this study is when the 
two points are on opposite boundaries, for instance r 2 
in region 1 and ri in region 3. In this case we found 
that all Green functions are of order 0{m^ ) 1 ). Then, the 
correlation functions given by 



(4.57) 



Takin g into acc ount all contributions from Eqs. (4.40), 
( 4.48 ) and ( 4.4£ ) to the total charge correlation function 
one finds 



/Li 
J S tota} (r 1) r 2 )dx 1 dx 2 



1 



2ttV 



will be finite of order O(l). But since there exists a 
surface charge density at the boundaries, the interesting 
quantity here is the surface charge correlation function 
{o~ s o~ s i) = dp, which will vanish in the limit S — > 0. 
We have the interesting result: the surface charge densi- 
ties at opposite boundaries are completely uncorrelated. 
To conclude this section let us return to the sum 



sinh 2 mW + 2a cosh mW sinh mW + a 2 sinh mW 



2ttV 



[cosh mW + a sinh mW] 
1 

1 j 

[cosh mW + a sinh mW] 



(4.58) 



rule (4.24) for the correlations functions along the bound- 
ary for W — > oo. If one takes into account all contri- 
butions (both particles in the fluid, one particle in the 
boundary and one in the fluid and both particles in the 
boundary) to the charge density correlation function, this 
function should read 



S , totai(ri,r 2 ) = <S(ri,r 2 ) 



(4.54a) 
(4.54b) 
(4.54c) 

<5(xi)<5( a ; 2 )( ( 7( a : 1 )a( a : 2 )) T (4.54d) 



S(xx) (a(xi)p(r 2 )y 
5{x 2 ) (cr(x 2 )p(ri)) J 



The sum rule (4.56) is not verified. The discrepancy how- 
ever is exponentially small when W — > oo. This situa- 
tion also occurs in the case a — which was studied in 
Ref. jy]. The two-component plasma is no longer in its 
conducting phase at T — 2 when confined in a slab. This 
is a special property of the two-component plasma at it 
is not related to the short-range attractive pot entia l near 
the boundaries. For this reason the sum rule ( 4.56| ) is no 
longer valid. Actually for a Coulomb syste m in a dielec- 
tric phase, the rhs. of the sum rule ( [4.56 ) should read 



(— l//37r 2 y 2 )(l — e^ 1 ) with e the effective static dielectric 
constant of the system. Then, in our case the system has 
a effective dielectric constant given by 



where the structure function S(rx, r 2 ) in Eq. ( 4.54a ) con- 
tains only the contributions for particles in the fluid, the 
correlations {po~) T in Eqs. ( 4.54b| ) and ( 4.54c ) contain 
the contribution when one point is on the boundary and 
the other in the film, and finally the correlation (o~o~) T 



(cosh mW + a sinh mWY 



(4.59) 



This phenomenon is particular to the two-dimensional 
two-component plasma and probably should not apply 
to real three-dimensional soap films. 
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V. CONCLUSION 

We have studied a toy model for electrolytic soap films. 
Although this model is very simple and gives only quali- 
tative information for real soap films it is very interesting 
since it is a solvable model. Our study of the disjoining 
pressure shows that the charging of the boundaries is 
responsible of the stability of the film. For strong adhe- 
sivity a > 1 the film is stable while for weak adhesivity 
a < 1 large films are not stable. For < a < 1 large 
films collapse to non-zero width small films. This could 
be the equivalent of a collapse to a Common Black Films 
for our two-dimensional model. For a — unstable large 
films collapse to a film of zero width which could be the 
equivalent of a Newton Black Film. We can conclude 
that the Coulomb interaction plays indeed an important 
role in the stability of large films. This is also the case 
in the one-dimensional model presented in Ref. [|| . Then 
it is natural to expect that for real three dimensional 
films the Coulomb interaction also plays an important 
role in their stability. Of course in real films there are 
other important interactions that certainly play a role in 
the stability of the film and in particular in the stability 
and structure of black films (both Common Black Films 
and Newton Black Films) that have not been taken into 
account in our simplified model. 

We also studied the density profiles and correlation 
functions for this two-dimensional model. The density 
profile near a boundary shows a classical double layered 



structure. A fraction of the anions (soap molecules) are 
stucked in the boundary creating a first layer of negative 
surface charge density. The ions in the fluid create the 
second layer of positive charge and thickness given by the 
screening length which screens the first layer. 

The correlation functions exhibit the usual behavior. 
In the x-direction across the film they decay exponen- 
tially with the characteristic screening length. In the 
y-dircction parallel to the boundary they decay alge- 
braically as 1/y 2 . The total charge correlation function 
(taking into account all contributions of particles in the 
fluid and in the boundary) obeys the usual sum rule for 
Coulomb fluids near a plane wall when W — > oo. For W 
finite, the two-component plasma at T — 2 is no longer a 
conductor and therefore fails to satisfy a sum rule for cor- 
relations along the boundaries. We also found an inter- 
esting new fact: the surface charge densities on opposite 
boundaries are completely uncorrelated. 
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